Two-dimensional Graphene with Structural Defects: 
Elastic Mean Free Path, Minimum Conductivity and Anderson Transition 
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Quantum transport properties of disordered graphene with structural defects (Stone- Wales and 
divacancies) are investigated using a realistic 7r — n* tight-binding model elaborated from ab initio 
calculations. Mean free paths and semiclassical conductivities are then computed as a function 
of the nature and density of defects (using an order-N real-space Kubo-Greenwood method). By 
increasing of the defect density, the decay of the semiclassical conductivities is predicted to saturate 
to a minimum value of 4e 2 /irh over a large range (plateau) of carrier density (> 0.5 10 14 cm -2 ). 
Additionally, strong contributions of quantum interferences suggest that the Anderson localization 
regime could be experimentally measurable for a defect density as low as 1%. 

PACS numbers: 73.23.-b, 72.15Rn, 73.43.Qt 



Clean graphene exhibits unique transport properties 
at the Dirac point. The density of charge carriers van- 
ishes but the conductivity remains finite in the order of 
a few e 2 /h [l], Q. This minimum of conductivity ob- 
served in the ballistic regime and related to contacts ef- 
fects is also present in the diffusive regime where disorder 
plays an important role. In conventional two-dimensional 
disordered metals, it is well established that such a low 
conductivity leads to an Anderson-type insulator at low 
temperatures 0, 01- In the Anderson regime, the elec- 
tronic states are spatially localized over a characteristic 
length scale, called localization length £ which drives 
the exponential suppression of conductivity with system 
length. The theory of Anderson localization initially de- 
veloped for electrons has been demonstrated to be ubiq- 
uitous in Physics, applicable to other types of particles, 
such as photons and atoms 0]. Surprisingly however, 
to date the observation of such metal-insulator transi- 
tion remains elusive in graphene, even in low-mobility 
devices. One of the reasons stems from the peculiar 
nature of localization effects in graphene whose charge 
carriers behave as masslcss Dirac fcrmions with chiral- 
ity degree of freedom at the origin of a sign reversal of 
the quantum correction to the semiclassical conductiv- 
ity (weak antilocalization When intervalley scatter- 
ing strongly predominates, ordinary weak localization is 
however expected to drive the system to an Anderson 
insulator [fl @] , but those theoretical results restrict to 
simplified disorder models, preserving the sp 2 symmetry. 

On the other hand, controlled defect engineering in 
sp 2 -carbon-based materials has become a topic of great 
excitment 0]. Indeed, the electronic (and transport) 



properties of carbon nanotubes [§| and graphene-based 
materials 13, U| can be considerably enriched by chem- 
ical modifications, including molecular doping and func- 
tionalization. Strong modifications of sp 2 -bondcd carbon 
materials by incorporation of sp 3 defects is also a suitable 
route for making graphene more sensitive to localization 
phenomena. Recently an attempt to turn graphene into 
a true band insulator was endeavored by transforming 



sp 2 bonds into sp 3 by hydrogenation 12|. Similar re 



suits have been experimentally debated in fluorinated 13 1 
or ozone treated graphene flakes 14 1. Using ion irra- 
diation, specific types of structural defects (vacancies) 
can be introduced in sp 2 -based carbon nanostructures. 
Convincing room-temperature signatures of an Anderson 
regime in Ar + irradiated carbon nanotubes have been re- 
ported (lij . In contrast, the conductivity of irradiated 
two-dimensional graphene saturates at the Dirac point 
above e 2 /h even down to cryogenic temperatures flij ]. 
suggesting a strong robustness of defective graphene. 

In this Letter, the electronic and transport proper- 
ties of disordered graphene are theoretically explored 
by introducing structural defects (Stone- Wales and di- 
vacancies) randomly distributed in the honeycomb lat- 
tice. The presence of these structural defects triggers 
resonant impurity levels, which broaden and generate im- 
purity bands as their density is increased. Ab initio cal- 
culations are performed to accurately describe the local 
energetics around the defects, and a tight-binding (TB) 
model (7T-7T* ) is elaborated from the corresponding band 
structures. Additionally, an efficient order-N real-space 
Kubo-Greenwood transport method is used to follow the 
transition from the diffusive to the insulating regime. 
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FIG. 1: (color online). Top panel : Schematic of various defects: Stone- Wales (a), 585 (b) and 555-777 (c) divacancies. Bottom 
panel: Corresponding electronic band structures computed using the SIESTA package (red lines) and parametrized using a 3 rd 
nearest-neighbors TB tt-tt* model (black dotted lines) along high-symmetry lines in the Brillouin zone (inset). Fermi energy is 
set at zero and the Dirac energy is indicated with a horizontal blue dashed line. 



In the present study, three types of non-magnetic 
structural defects are considered. The first one is the 
so-called Stone- Wales (SW) defect which consists of two 
heptagons connected with two pentagons (FigUJa). Al- 
though its presence has already been experimentally re- 
ported in graphene (Io| . its influence on the transport 
properties remains unknown. Ma et al. have reported 
in a theoretical studyQ that a SW defect could produce 
a slight out-of-plane deformation. We consider here an 
in-plane geometry as a first approach. The two other de- 
fects consist in various reconstructions of the divacancy : 
one being associated to the formation of 2 pentagons and 
1 octagon (585 - FigJTJb), while the other is composed of 
3 pentagons and 3 heptagons (555-777 - Fig[TJc). Our ab 
initio calculations predict the 555-777 reconstruction to 
be more stable than the 585 one by ~ 0.9 eV. in agree- 
ment with previous theoretical predictions jl7| . Finally, 
divacancies are known to be more stable than two iso- 
lated monovacancies whose migration energy barrier is 
rather low [81. fl7j|. In contrast to the monovacancies, the 
divacancies do not require a spin-dependent electronic 
structure treatment [181 ]. 

To reduce the computational cost of the transport 
calculations, the electronic band structures (FigJIJ first 
computed with the SIESTA code [H, H(j are then re- 
produced using a parametrized 3 rd nearest-neighbors ir- 



7r* TB Hamiltonian model [22|, |23[ . The agreement be- 



tween this TB model and the ab initio calculation is 
quite satisfactory, especially in the transport energy re- 
gion [— l,l]eV around the Dirac point. The observed 
differences at higher energies are basically related to the 
inability of the 7r-7r* TB model to reproduce the con- 
duction bands of graphene along K-M branch. The SW 
defect does not display any doping character since both 
Fermi and Dirac energies are aligned (Fig[TJa). On the 
contrary, both divacancies are found to act as acceptors 



since Fermi energy is below the Dirac point (Fig[TJb-c). 

Based on the present TB model, the density of states 
(DOS) are calculated using the recursion method, which 
is very efficient for large disordered systems [24| . In FigJ5] 
(inset), the total DOS corresponding to a defect density 
(rid) of 1% are displayed for the three different defects. 
A large difference in the position of the defect-induced 
resonant impurity levels can be observed. For instance, 
the quasi-bound states localized around the SW defect 
yields a bump in the DOS around 0.35eV above the Dirac 
energy (Figj2]- solid line). The 585 defect exhibits a sim- 
ilar DOS fingerprint but with an opposite behaviour to 
the case of SW defect. Indeed, for this defect the corre- 
sponding peak appears around — 0.35eV (FigJ5]- dashed 
line). Finally, in presence of 555-777 defects, two differ- 
ent DOS impurity-peaks are reported : a wider impurity 
band in the hole region (— 0.8eV) and a sharp peak on the 
electron side (0.6eV) (Fig(5]- dashed-dotted line). In all 
cases, the intensity and the width of the defect-induced 
peaks are found to increase with rid- 

Transport properties of large and disordered graphene 
systems are then calculated using an efficient order- 
N Kubo-Greenwood method 24|. This powerful tech- 



nique gives a direct access to the elastic mean free paths 
at energy E, extracted from the wavepacket dynamics. 
The latter is characterized by the time-dependent dif- 
fusivity coefficient D(E,t) = AR 2 (E,t)/t with AR 2 = 
AX 2 + AY 2 and AX 2 {E 1 t) = Tr[<J(E - H)|X(t) - 
X(0)| 2 ]/Tr[5(E - H)]. Tr is the trace over p z orbitals 
and Tr[5(E - H)]/S = p(E) is the total DOS (per 
unit of surface). The two position operators X(t) and 
Y(t) are expressed in the Heisenberg representation 
(X(t) = W(t)X(0)U(t)) and the time evolution oper- 
ator U(t) = 11^=0 exp(iHAt/h) with At the chosen time 
step, is computed with a Chebyshev polynomial expan- 
sion method 2J|. Calculations are performed for several 
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FIG. 2: (color online). Elastic mean free paths (main frame) 
and total DOS (inset) for disordered graphene with a defect 
density of 1% of SW (solid line), 585 (dashed line) or 555-777 
(dotted dashed line) divacancies. All curves have been aligned 
to make coincide the minimum of DOS with the zero energy. 
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FIG. 3: (color online). Semiclassical conductivities in dis- 
ordered graphene versus carrier densities and various defect 
densities (from 0.1% to 1%) of Stone-Wales (left panel), 585 
(right panel main frame) and 555-777 (right panel inset). 



initial random phase wavepackets, and for total elapsed 
time t « 1.5ps. The typical size of the simulated system 
is ~ 0.074/im 2 (containing 2.8 x 10 6 carbon atoms), large 
enough to avoid finite size effects. In the Kubo formalism, 
the different transport regimes can be inferred from the 
behaviour of D(E,t). The wavepacket velocity v(E) can 
be extracted from the short time behavior of the diffu- 
sivity, D(E, t) ~ v 2 (E)t, while the elastic mean free path 
£ e {E) is estimated from the maximum of the diffusivity, 
fl mal (£) = 2v(E)£ e (E). Finally the Kubo semiclassical 
conductivity reads a sc {E) = \e 2 p(E)D max (E) . 

The mean free paths are calculated for defect densities 
varying from rid = 0.1% to 1%. At a given energy, £ e is 
predicted to be inversely proportional to rid, as expected 
from the Fermi golden rule. £ e also displays a strong en- 
ergy dependence, with dips around the Dirac point where 
it reaches values as low as few nanometers. These dips 
are associated with the resonant impurity states, which 
are specific to each defect. £ e (E) can change by one or- 
der of magnitude around the Dirac point depending on 
the nature of the defect and the energy of charge carri- 
ers, as shown on Fig(2] (main frame) for a defect density 
rid = 1%. For instance, l e is estimated to be ~ 5nm 
at the Dirac point in presence of 585 defects, whereas 
£ e ~ 50nm for the 555-777 defects because the two as- 
sociated energy resonances (—0.8 and 0.6eV) are farther 
from the Dirac point. Consequently, this defect nature 
dependence will strongly impact both semiclassical and 
quantum transport regimes, as discussed below. 

Figure [3] illustrates the semiclassical conductivities for 
disordered graphene with the three possible defects eval- 
uated as cr sc (n) = ie 2 p(n)v(n)£ e (n) , where n(E) is the 
carrier density |25| . In addition to the strong dependence 
of a sc on the carrier energy, a saturation of the conductiv- 
ity decay with rid is observed, reaching the minimal value 
of <7™ ln = 4e 2 jixh. As rid increases, a large plateau of 



minimum conductivity develops close to the Dirac point, 
similarly to the case of simpler point defects [26|. The 
origin of such a minimum of conductivity has recently 
been strongly debated in the literature |27| ■ Our simula- 
tions suggest that such a value could be associated to a 
graphene system containing structural defects, in absence 
of quantum interferences (QI). One further observes that 
for each type of defects, the plateaus of minimal con- 
ductivity extend over the energy window encompassing 
the corresponding resonant impurity bands. Indeed, for 
the SW case, this plateau is primarily located at ener- 
gies above the Dirac point, in contrast with the 585 case, 
where the minimum conductivity is observed below the 
Dirac point. The case of 555-777 is more complicated 
due to the presence of two impurity resonant energy win- 
dows. In the present study, a rigid defect model has been 
considered, ignoring the screening effects which could oc- 
cur at high carrier densities. However, we have carefully 
checked that increasing the carrier density in the ab initio 
calculations does not significantly modify the bandstruc- 
ture, supporting this rigid defect model approximation. 

One noteworthy observation is that the predicted short 
mean free paths close to the Dirac point favors strong 
contributions of QI, as long as the transport regime re- 
mains quantum coherent. In the present Kubo formal- 
ism, these QI contributions can be evidenced in the ratio 
D(t) I -Dmax which departs from unity at long times in the 
presence of QI. Fig0] (main panel) shows several typical 
behaviors of D(t) / D max at selected energies and for the 
case of rid = 1% of SW defects. A more global picture 
(over a larger part of the energy spectrum) is also given 
in the inset. At E = 0.5eV (close enough to the res- 
onance energy associated to the SW, see Fig. 2 (inset)), 
D(t) I 'Dmax exhibits a fast decay consistent with the es- 
timated short localization length (£ ~ lOnm). Indeed, 
following the scaling theory of localization, £ can be de- 
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FIG. 4: (color online). Time dependent diffusion coefficient 
D(E, t) normalized to D max = 2v£ e at three selected energies 
and for a 1% density of SW defects. Inset: Enlarged view of 
D(f)/_D max over the energy spectrum. 



tcrmined once the semiclassical transport length scales 
arc known, as £(E) = £ e (E) cxp(7r/icr sc (£')/2e 2 ) Q. Sim- 
ilar values are found at the resonant energies for both 1% 
of 585 defects (£ - lOnm) and for 1% of 555-777 defects 
(£ ~ 25nm). Actually, in the energy windows where a sc 
saturates to cr™ ln , short localization lengths (£ < lOOnm) 
are obtained for all defect cases. Much longer £ are 
predicted for energies outside the minimum conductivity 
plateaus, as illustrated in Fig. 4 for E = 1.25cV where the 
contribution of QI can even become vanishingly small. 

Specific conductance fingerprints will be thus obtained 
depending on the nature and density of defects. Although 
real samples of defective (irradiated) graphenc are likely 
to encompass a mixture of those different defects, the sat- 
uration of the semiclassical conductivities and the typical 
localization lengths (in the range of 30nm for 1% de- 
fects) close to the Dirac point will be a robust common 
feature to all possible cases. Considering the reported 
weak electron-phonon coupling 28j , such insulating state 
should be observable even at room temperatures (using 
suspended graphenc preferentially). 
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